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. Abstract. The definition of the standard differential operator is extended from 

integer steps to arbitrary stepsize. The classical, nonrelativistic Hamiltonian is 
quantized, using these new continuous operators. The resulting Schroedinger type 
equation generates free particle solutions, which are confined in space. The angular 
momentum eigenvalues are calculated algebraically. It is shown, that the charmonium 
spectrum may be classified by the derived angular momentum eigenvalues for 
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1. Introduction of the continuous differential operator 

> : 

Since Newton PQ and Leibniz [2] introduced the concept of infinitesimal calculus, 
differentiating a function with respect to the variable x l is a standard technique applied 
in all branches of physics. The differential operator di, 

transforms like a vector, its contraction yields the Laplace-operator 

A = d% (2) 

which is a second order derivative operator, the essential contribution to establish a 
wave equation, which is the starting point to describe several kinds of wave phenomena. 

Until now, the differential operator has only been used in integer steps. We want 
to extend the idea of different at ion to arbitrary, not necessarily integer steps. A natural 
generalization is to search for an operator D by setting 

D m = d n (3) 

where m, n are integers. Formally, this is solved by extracting the m-th root 

D = d n/m m,neN (4) 
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or, even more general, we will introduce a continuous differential operator by 

D = d p P^R+ (5) 

with p beeing a positive, real number. 

In order to derive a specific representation for the operator D, we introduce a very 
simple set of basis functions: 

f(x) =x v veR + (6) 

Differentiation of these functions is well defined for the integer case. It is just a 
multiplication with an appropriate factor and a linear reduction of the exponent. 

This procedure can easily be extended to the continuous case. 

The solution is inspired by a similar problem, which already has been solved in 
mathematical history: the extension of the interpretation of the faculty function n\ to 
real numbers, which resulted in the famous Euler Gamma function Y{x + 1), which is 
widely used in many different of analysis. 

Thus, the definition for the continuous derivative is: 



Dx u = d p x v 



for v = and p > 

V(u + 1) (7) 
; x u - p forz/,p>0, u-p>0, x>0 



[T(u+l-p) 

The set of functions {x v } is sufficient to span a Hilbert space. 
Using the above definition, we are able to quantize the Hamiltonian of a classical 
particle and solve the corresponding Schroedinger equation. 

2. Quantization of the classical Hamiltonian and free particle solutions 

We define the following set of conjugated operators 

P^ = {P ,P t } ={ihd t ,-ihD i } (8) 
X II = {X ,X t } = {t,x P /T(p+l)} (9) 

which satisfy the following commutator relations 

[x,, x 3 ] = o, [Pi, p 3 ] = o, [x,, p 3 \ = -ihSij (io) 

With these operators, the classical, non relativistic Hamilton function H c 

H.= ^ + V(«0 (ID 
is quantized. This yields the Hamiltonian H 

H = ^—D i D i + V(X*) (12) 

Thus, a time dependent Schroedinger type equation for continuous differential operators 
results 

Hm = (—DWi + VtX*)) * = ihdt* (13) 
\ 2m / 




Figure 1. Free particle solutions for the continuous differential operator Schroedinger 
type equation, for different values of p, p — 1, 0.9, 0.8, 0.7, 0.6. Solid lines are cos(p, x), 
dashed lines sin(p, x). For p = 1 solutions reduce to the standard cos(x) and sin(x) 
functions. For p < 1 these functions are increasingly located at x = 



For p = 1 this reduces to the classical Schroedinger equation. 

We will now present the free particle solutions for this equation. We can do this, 
since for V(X l ) = the commutator [P^, H] vanishes and consequently, energy and 
momentum are conserved. 

We extend the definition for the standard series expansion for the sine and cosine 
function 

oo 

sm(p, x) = J2 (-l) n x {2n+1)p /T((2n + l)p + 1) (14) 

n=0 

oo 

cos(p,x) = ^(-l)V" p /r(2np+ 1) xeR+ (15) 

n=0 

With this definition, the following relations hold 

D sin(p, kx) = k p cos(p, kx) (16) 
D cos(p, kx) = — k p sin(p, kx) k,x G R + (17) 

It follows from these relations, that this functions are the eigenfunctions of the free 
Schroedinger type equation. In the stationary case we get the energy relation 

E = —k 2p (18) 

In figure 1 the functions sin(p, x) and cosQo, x) are plotted for different values of 
p. While for p = 1, these functions reduce to the known cos(x) and sin(a;), which 
are spread over the whole x-region, for p < 1 these functions become more and more 
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located at x = and oscillations are damped, a behaviour, which we know e.g. from 
the Airy-functions. 

Let us make this observation more plausible: for any free particle solution w(p) 
with p = 2/m and m odd, which obeys the differential equation 

D 2 w = -k 2p w (19) 

differentiating m — 2 times we obtain: 

D m -\D 2 w) = D m - 2 (-k 2p w) (20) 
D m w = -k 2p D m - 2 w (21) 

d 2 w = ±k 2p k^ m -^ p (—^—— + R(p, x))w, lim R(p, x) = Q (22) 

r(l + p) 

= ±k 2 f ^- ) w + o(x^) (23) 

This corresponds to an additional, impulse dependent potential with a leading linear 
term for similar solutions of the classical Schroedinger equation (p = 1). 

Thus, as the main result of our derivations so far, free particle solutions for p ^ 1 
are localized in space. Consequently, if there are particles, which show such an behaviour 
in nature, a Schroedinger equation with continuous derivatives could be a useful tool 
for a description. 

In order to obtain more properties of the continuous differential operator 
Schroedinger equation, we will now calculate the eigenvalues of the angular momentum 
operator. 

3. Classification of angular momentum eigenstates 

We define the generators of infinitesimal rotations in the i, j-plane (i,j = 1, ...,37V), TV 
is the number of particles): 

= XiPj - XjP { (24) 

Since [Lij,H] = 0, angular momentum is conserved. Commutator relations for are 
isomorph to a SO(3N) algebra: 

[Lij, Zj mn ] Hx(5i m L j n -\- d~j n Li m d~i n Lj m 5j m Li ri ^ (26) 

Consequently, we can proceed in a standard way j^j, by defining the Casimir-Operators 

A^EOM 2 , k = 2, ...,37V (27) 

which indeed fulfill relations [K\ Nl Lif\ = and successively [A|, = 0. Their explicit 
form is given by 

A 2 k = -tf^Di + X^iDWj - SlXWj - X'X'DiDj) (28) 
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Now we introduce a generalization of the homogenous Euler operator for continuous 
differential operators 

J k (p)=X l D l (29) 
= (xY/T(p+l)D t (30) 

With the generalized Euler operator the Casimir-operators are: 

A* = -tf{J k e + X^DW, - U k e - J%J k e - I)) (31) 

Now we define a Hilbert space H p of all homogenous polynoms of degree IkP, which 
satisfy the Laplace equation D l Dif = 0: 

H p = {/ : f(\x) = X lkP f(x); DWJ = 0} (32) 

On this Hilbert space, the generalized Euler operator J k (p) is diagonal and has the 
eigenvalues 

' for n = 



h(p,n) 



^ + for, = 1,2,3,... (33) 



{ r(p + i)r((n-i) P + i) 

We want to emphasize, that these eigenvalues are different from the degree of homogenity 
in the general case p ^ 1, or, in other words: only in the case of p = 1 homogenity degree 
n of the polynoms considered coincides with the eigenvalues of J k (p = l,n). 

Once the eigenvalues of the generalized Euler operator are known, the eigenvalues 
of the Casimir-operators A 2 , are known, too: 

A 2 f = M 2 (p,n)f (34) 

A 2 J = h%(p,n)(l k (p,n) + k — 2)f (35) 

with 

k > 4-i > ... > k > (36) 

For the case of only one particle (N = 1), we can introduce the quantum numbers 
j and m, which now denote the j-th or m-th eigenvalue of the Euler operator. The 
eigenfunctions are fully determined by these two quantum numbers / =| jm > 
With the definitions L z = L 12 and J 2 = L\ 2 + Lj 3 + L| 3 it follows 

L z | jm > — hkipjin) | jm > m — 0, +1, +2, +j (37) 

J 2 \jm> =h%(p,j)(l 3 (p,j) + l)\jm> j = 0,+l,+2,... (38) 

Please note the fact, that in the case p ^ 1 L z has no negative eigenvalues. 

In table 1 the first seven eigenvalues of L z and J 2 for a single particle are listed for 
p = 1 and p = 2/3. For n > 1 and p < 1 the eigenvalues of the Euler operator are not 
equally spaced any more, instead stepsize is strongly reduced. Since the Euler operator 
eigenvalues contribute quadratically into the definition J 2 , the energy of higher total 
angular momenta is reduced increasingly. 
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Table 1. Eigenvalues in units of h for the angular momentum states of a single particle. 
n is a counter for the eigenvalues of the Euler operator, eigenvalues for L z and J 2 are 



given for p 


= 1 and 


p = 2/3. 






n L z (p = 


l,m) 


L z (p = 2/3, m) 




J 2 (P = 2/3, j) 
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1 


2 
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2 


1.460 998 49 


6 


3.595 515 06 


3 


3 


1.860 735 02 


12 


5.323 069 84 


4 


4 


2.222 222 22 


20 


7.160 493 83 


5 


5 


2.556 747 35 


30 


9.093 704 37 


G 


6 


2.870 848 32 


42 


11.112 618 39 



We have derived the full spectrum of the angular momentum operator for the 
continuous differential operator Schroedinger type wave equation by use of standard 
algebraic methods. 

We expect additional information about the properties of this wave equation, if we 
consider its linearized pendant. We present some results in the next section. 

4. Results for the linearization of a second order differential equation 

Linearization of a second order wave equation was first considered by Diracjl]. We 
propose a generalization of his approach, starting with a linear version of the wave 
equation, which is iterated m times, roughly: 

( 7 D) m = d 2 (39) 

This is formally solved by taking the m-th root: 

1 D = (l) (1/m) 9 2/m (40) 

A new property is revealed via linearization: a phase factor, which may not be neglected. 
If m is an integer, m different phase factors exist, which have to be incorporated within 
the definition of 7. Consequently, 7 may be interpreted as a matrix with at least mxm 
dimensions. 

The 7 matrices obey an extended Clifford algebra (P means all permutations): 

m 

EIl7i = (41) 

P i=l 

This equation is fulfilled by a set of hermitean, traceless mxm matrices, which built 
at least a subspace of SU(m). 

Finally, we can distinguish two distinct parts of linarization. If spacelike 
components are linearized (D = Di), an additional SU(m) symmetry is introduced, 
which has been demonstrated for the case m = 2 by Levy-Leblond by linearizing 
the ordinary Schroedinger equation to get an additional SU(2) symmetry jHj. For a 
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single linearized timelike component, we obtain m different energy values, which are 
interpreted as m different particles. 

From a mathematical point of view, a continuous differential operator may be 
interesting for any value of p, a physicist may first concentrate on p = 2/m, because 
there is a simple interpretation of the physical properties described by such a wave 
function. 

We want to emphasize, that p = 2/3 seems to be the most promising candidate for 
a physical interpretation: 

It matches a triple iteration of a linear wave equation, which shows a direct 
connection to SU(3) and provides a description for up to 3 particles simultaneously. 
Furthermore, as we have already shown, particles are automatically confined in space. 

Summarizing all these facts, we assume, that the continuous differential operator 
Schroedinger type wave equation with p = 2/3 is an appropriate candidate for a non 
relativistic description of particles with quark-like properties. 

5. Interpretation of the Charmonium spectrum 

In the previous sections we have introduced the idea of continuous differential operators 
and discussed some properties of the resulting non relativistic Schroedinger equation 
and its linearized pendant. 

We have developed a new theoretical concept, which fulfills at least the following 
three demands: First, available experimental data will be reproduced with a reasonable 
accuracy. Second, it will give new insights on underlying symmetries and properties of 
the objects under consideration. Third, we will make predictions, which can be proven 
by experiment. 

None of our results, presented so far, does require any information from QCD or a 
similar theory. All our statements could have been made in the 1930ies already, even 
though they would have been highly speculative. Today, we are in the comfortable 
position, that there are enough experimental data, our predictions can be compared 
with. 

A promising candidate is the charmonium spectrum jS]. 

In the upper part of figure 2 we have displayed all experimentally observed 
charmonium-states with experimental masses, which are normally compared with results 
from a potential model, which tries to simulate confinement and attraction by fitting a 
model potential |Zj,|H]- 

We will assume, that this spectrum is a single particle spectrum for a particle, 
whose properties are described by the continuous differential operator Schroedinger 
equation ((TBI with p = 2/3. We suppose, the system is rotating in a minimally coupled 
field, which causes a constant magnetic field B. This leads to the following Hamiltonian 
H or mass formula 

H | jm>= m c 2 + a J 2 + BL Z \ jm > (42) 
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Figure 2. Charmonium spectrum. All observed particles are given with their name, 
experimental mass, the proposed SO (3) conform quantum numbers and calculated 
mass according to mass formulas l|42ll . (|47|l . States are labeled with the derived SO (3) 
quantum numbers, with j and m are the j-th and m-th eigenvalue of the generalized 
Euler operator. For j — 3, the influence of an additional constant correction term At 
is illustrated. 



where moc 2 , a and B will be adjusted to the experimental data. The eigenf unctions are 
modified spherical harmonics with D l Di \ jm >= 0, the eigenvalues for J 2 and L z are 
given by (}3Tj) . (j38|) and are listed in table 1. 

In figure 2 we have sketched the resulting level scheme and the quantum numbers 
< jm > for every particle. 

We have to prove, that this leads to correct results. 

The first crucial test is the verification of the correct value of the non trivial m = 2 
quantum number which corresponds to the n = 2 eigenvalue of the generalized Euler 
operator. For the set of x-particles, we obtain: 

T < ■ 9 9^ X(3556) < 22 > -x(3415) < 20 > 

L z (j = 2, m = 2) exp = x(35n)<21> _ x(3415)<20> = 1-468 (43) 

This is an exact match with the theoretical value within the experimental errors, which 
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are about ±2[MeV}. 

Next we will determine the constants m c 2 and a. We choose the two lowest 
experimental states of the standard charmonium spectrum, ?7 C (2980) < 10 > and 
x(3415) < 20 >. This yields a set of equations 

m c 2 + 2a = ry c (2980) < 10 > (44) 

m c 2 + 3.59551506a = x(3415) < 20 > (45) 

which determine m c 2 = 2434.72[MeV] and a = 212M[MeV] 

Our level scheme predicts a particle with quantum numbers < 00 >, which is 
beyond the scope of charmonium potential models. According to our mass formula, it 
has a predicted mass of m c 2 = 2434. 72 [Me V). 

This is the second crucial test, since this is a low lying state, it should already have 
been observed. Indeed, there exists an appropriate candidate, the Ci(E c )(2430) < 00 > 
particle [H|,[H|. 

This result supports our theory, since it is a second occurance of the m = 2 
eigenvalue of the generalized Euler operator. 

In addition, it supports the idea, that the energy levels of the spectrum are labeled 
correctly. 

Finally, the minimal difference of only 4[MeV], which is mainly due to the 
uncertainty of the ?y c (2980) < 10 > experimental mass, between predicted and 
experimental mass of the Ci(E c )(2430) < 00 > particle indicates, that the assumed 
SO (3) symmetry is fulfilled exactly. 

As a last step, we need to determine the value of B. From the experimental data, 
we observe a slight dependence on j, which we will ignore for sake of simplicity. Instead, 
we do a least square fit for B(j = 1) and B(j = 2) and obtain B = 107[MeV]. 

Next step is to find higher eigenvalues of the generalized Euler operator. As a 
matter of fact, all candidates are positioned above 3.73[GeV], the threshold for DD- 
meson production. From the point of view of our theory, the particle properties change 
and we should actually make a new fit for this particle with Hp mass formula and new 
parameters. 

H D | jm >= m 0D c 2 + arjJ 2 + B D L Z \ jm > (46) 

For B])we obtain: 

B D = ^(4040) < 31 > -^(3770) < 30 >= 270[MeV] (47) 
B D = 2.5B (48) 

From a pessimistic point of view, this is just one more parameter, from an optimistic 
point of view we could associate a correlation of B and R, the ratio of hadronic and 
mesonic cross sections, which shows a similar behaviour |l()j. 
We first prove the position of the ^(4160) < 32 > state 

*(4160) < 32 > -W3770) < 30 > AlAA . inS 

LJj = 3, m = 2) exv = -) = 1.444 49 

U ' ' p (4040) < 31 > -tf (3770) < 30 > V ' 
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This is in exact agreement with the theoretical value within the experimental errors, 
which are about ±10 [Mel 7 ] and the third observation of the n = 2 eigenvalue of the 
generalized Euler operator. This is an indication, that the proposed level scheme is still 
valid in the above region. As a consequence, we expect the < 33 > state at 



This state is missing in the experimental spectrum, so this is the second predicted 
particle or resonance, respectively in this area. 

Consequently, we have only one experimental candidate for for j = 3 and j — 4 
respectively. With the assumption op = a and mo£>c 2 = tuqC 2 we obtain the theoretical 
values 



For j = 3 the calculated mass differs by 116 [MeV] from the experimental value. In 
other words, the j = 3 eigenvalue of the generalized Euler operator may be deduced 
from experiment as 1.77 instead of 1.86 from the theory. This is a deviation of about 
5%, which is close, but compared with the previous results not precise. 

On the other hand, the theoretical #(4410) < 40 > mass matches exactly with the 
experimental value within the experimental errors. 

This indicates, that the particles for j = 3, observed in experiment, carry an 
additional property, which reduces the mass by the amount of e.g. a pion. Of course, 
if we add an additional (Ar)5j 3 term to the proposed mass formula, we can shift these 
levels by the necessary amount. Figure 2 shows the theoretical masses for j = 3 with 
such an correction term switched off and on. 

Summarizing these results, the charmonium spectrum reveals an underlying SO (3) 
symmetry, which agrees with the predictions of our theory in the case of p = 2/3. 
The eigenvalues of the generalized Euler operator conform within experimental errors 
with experimental data for n — 0, 1, 2, 4, for n = 3 there is a deviation of 5%. Two new 
particles have been predicted. With only 4 parameters (including At) the massformulas 
(I42|) . (|47jl reproduce the charmonium spectrum within a range of 2GeV. 

6. Conclusion 

We have defined a continuous differential operator for arbitrary stepsize p. A 
Schroedinger type wave equation, derived by quantization of the classical Hamiltonian, 
bases on these operators, generates free particle solutions, which are confined to a certain 
region of space. The multiplets of the generalized angular momentum operator have been 
classified acoording to the SO (3) scheme, the spectrum of the Casimir-Operators has 
been calculated analytically. We have also shown, that for p = 2/3, corresponding to a 
triple iterated linearized wave function an inherent SU(3) symmetry is apparent. 



*(?) < 33 > = #(3770) < 30 > +1.86073502^ 
= 4272 [MeV] ± 15 [MeV] 



(50) 
(51) 



#(3770) < 30 > th = 3886[MeV] 
#(4410) < 40 > th = 4387[MeV] 



(52) 
(53) 
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From a detailed discussion of the charmonium spectrum we conclude, that the 
spectrum may be understood within the framework of our theory. 

Two new particles have been predicted, one of them already confirmed by 
experiment, the other being a resonance at 4272 [MeV] ± 15 [Mel 7 ]. 

The results indicate, that continous differential operator wave equations may play 
an important role for our understanding of phenomena in the area of particles with 
quark-like properties, e.g. confinement. 

Up to now, we have only derived a non relativistic wave equation. A next step will 
include relativistic effects too. 
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